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Abstract 

For identical bosons with a large scattering length, the dependence of the 3-body recombina- 
tion rate on the collision energy is determined in the zero-range limit by universal functions of 
a single scaling variable. There are six scaling functions for angular momentum zero and one 
scaling function for each higher partial wave. We calculate these universal functions by solving 
the Skorniakov~Ter-Martirosian equation. The results for the 3-body recombination as a func- 
tion of the collision energy are in good agreement with previous results from solving the 3-body 
Schrodinger equation for ^He atoms. The universal scaling functions can be used to calculate the 
3-body recombination rate at nonzero temperature. We obtain an excellent fit to the data from 
the Innsbruck group for ^^^Cs atoms with a large positive scattering length. 
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I. INTRODUCTION 



Three-body recombination is a 3-particle process in which two of the particles bind to 
form a molecule. This process is important in cold atom physics, because it is one of the 
most important loss mechanisms for trapped ultracold atoms. If the interactions of the 
atoms are known with sufficient accuracy, the 3-body recombination rate can be calculated 
by solving the 3-body Schrodinger equation numerically. However the interactions between 
complex atoms may not be known with sufficient accuracy. This is especially true for atoms 
whose S-wave scattering length a is large compared to the range of their interaction. In this 
case, the atoms have universal properties that depend on a but are otherwise insensitive 
to length scales set by the range. This makes possible a complementary approach to the 
few-body problem that involves expanding in powers of the range divided by a. The leading 
term in this expansion is called the zero-range limit or, alternatively, the scaling limit. In 
the 3-body sector, the approach to the few-body problem based on expanding around the 
zero-range limit was pioneered by Vitaly Efimov beginning around 1970. 

Efimov discovered that the zero-range limit of the 3-body problem for nonrelativistic 
particles with short-range interactions is quite remarkable. If a = icso, there are infinitely 
many 3-body bound states with an accumulation point at the 3-atom scattering threshold. 
These Efimov states or Efimov trimers have a geometric spectrum [1] : 



where is the binding wavenumber of the Efimov trimer labeled by and m is the mass 
of the particles. This geometric spectrum is a signature of a discrete scaling symmetry with 
discrete scaling factor e'^/*" [2]. In the case of identical bosons, Sq ~ 1.00624 and the discrete 
scahng factor is e'^^^° ^ 22.7. We will refer to the universal phenomena characterized by 
this discrete scaling symmetry as Efimov physics [3]. Efimov showed that discrete scale 
invariance is also relevant if a is large but finite [4, 5]. The negative values of a for which 
there is an Efimov trimer at the 3-atom scattering threshold are related by the discrete 
scaling symmetry. The positive values of a for which there is an Efimov trimer at the 
atom-dimer scattering threshold are also related by the discrete scaling symmetry. Another 
dramatic example of Efimov physics is the vanishing of the 3-body recombination rate at 
threshold at positive values of a that are related by the discrete scaling symmetry [6-8]. 

Efimov also showed that the discrete scaling symmetry governs the energy dependence of 
3-body scattering processes, such as atom-dimer elastic scattering and 3-body recombination 
[5]. It is convenient to measure the total energy E of the three atoms in their center-of-mass 
frame relative to the 3-atom scattering threshold. The discrete scaling symmetry implies 
that the scattering rates at energy in a system with parameters a and differ from those 
at energy (e~^'^/*°)"£^ in a system with parameters (e'^/*°)'*a and only by an overall change 
in the scale. The zero-range limit gives more detailed predictions that Efimov referred to 
as the Radial Law. The Radial Law expresses the dependence of the scattering rates on E 
in terms of universal functions of a scaling variable x = {a^mE jU^y/'^ . Once these scaling 
functions have been calculated, the scattering rates can be predicted for any system of 
identical bosons with a large scattering length. 

The universal results defined by the scaling limit are exact only in the limit of zero range. 
Corrections to the universal results associated with a nonzero range are suppressed by powers 
of £/|a|, where I is the natural low-energy length scale [2]. For a short-range potential, £ 
is simply the range of the potential. For atoms whose potential has a van der Waals tail 
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—Cq/t^^ t is the van der Waals length £vdw — {jnC^jTi^'Y^^ ■ The effects of the nonzero range 
can be calculated within a systematic expansion in the small parameter ^/|a| [9-13]. The 
next-to-leading order term in this expansion is proportional to r^/ a, where is the effective 
range [9, 10]. Surprisingly, the next-to-next-to-leading order term in this expansion is also 
determined only by the scattering length a, the Efimov parameter k*, and the effective range 
[13]. 

The natural low-energy length scale i allows us to differentiate between a shallow dimer 
and deep dimers. The shallow dimer is a universal feature of the two-body system with 
large positive scattering length and has binding energy ~ Ti^/{ma'^). The underlying 
interaction can also support deep dimers whose binding energies are comparable to or larger 
than h'^/(m£'^). The properties of these deep dimers are not universal. The alkali atoms 
used in most cold atom experiments form diatomic molecules with many deeply-bound 
energy levels. Efimov trimers can therefore decay into an atom and a deep dimer with large 
kinetic energies. These final states also provide inelastic atom-dimer scattering channels 
and additional 3-body recombination channels. If the atoms do not form deep dimers, the 
zero-range results for scattering rates are determined by two parameters: the scattering 
length a and the Efimov parameter k^. If the atoms do form deep dimers, the zero-range 
results for scattering rates are determined by three parameters: a, k*, and a parameter 77* 
that determines the widths of the Efimov trimers [14, 15]. Remarkably, the same universal 
functions that determine the energy dependence of scattering rates for the case 77^ = in 
which there are no deep dimers also determine their energy dependence for 77* > 0. 

The first experimental evidence for Efimov physics was presented by the Innsbruck group 
[16]. They carried out experiments with ultracold ^'^'^Cs atoms in the lowest hyperfine 
state, using a magnetic field to control the scattering length. They observed a resonant 
enhancement in the 3-body recombination rate at a ~ —850 Gq that can be attributed to 
an Efimov trimer near the 3-atom threshold. At the temperature 10 nK, the loss rate as a 
function of a can be fitted rather well by the universal formula for zero temperature derived 
in Ref. [15] with a width parameter 1]^ = 0.06(1). The Innsbruck group also observed a local 
minimum in the 3-body recombination rate near a ^ 210 ao that might be associated with 
an interference minimum of the 3-body recombination rate. The measurements were carried 
out at 200 nK, which is large enough that taking into account the nonzero temperature is 
essential. 

The 3-body recombination rate can be calculated at nonzero temperature by carrying 
out a thermal average of the 3-body recombination rate as a function of the collision en- 
ergy. There have been several previous calculations of the 3-body recombination rate for 
atoms with large scattering length at nonzero temperature using specific models. DTncao, 
Suno, and Esry used a simple 2-parameter potential with either one or two S-wave bound 
states [17, 18]. Massignan and Stoof used a 4-paramctcr scattering model for atoms near 
a Feshbach resonance [19]. There have been several previous calculations of the 3-body 
recombination rate for nonzero temperature that have exploited the universality of atoms 
with large scattering length. In the case a < 0, the 3-body recombination rate was calcu- 
lated by Jonsell using the adiabatic hyperspherical approximation [20] and by Yamashita, 
Frederico, and Tomio using a resonance approximation [21]. In the case a > 0, the 3-body 
recombination rate has been calculated using simplifying assumptions to neglect some of the 
universal scahng functions [22-24]. The calculations in Refs. [20-24] all involve uncontrolled 
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approximations, so they are not definitive zero-range predictions.^ 

In this paper, we present definitive zero-range predictions for the 3-body recombination 
rate at nonzero temperature for the case a > 0. We calculate the universal scaling functions 
that determine the hyperangular average of the 3-body recombination rate at collision ener- 
gies up to about 30 times the binding energy of the shallow dimcr. Our results are validated 
by comparing with a previous calculation of the 3-body recombination rate as a function 
of collision energy for ''He atoms. The universal scaling functions can be used to calculate 
the rate constants for 3-body recombination into the shallow dimer and into deep dimers at 
temperatures and densities where Boltzmann statistics are applicable. We apply the results 
to ^^^Cs atoms with large positive scattering length a and show that they give an excellent 
fit to the data from the Innsbruck group on the 3-body recombination rate as a function of 
a. 



II. THREE-BODY RECOMBINATION 

Three-body recombination is a 3-atom collision process in which two of the atoms bind 

to form a diatomic molecule or dimer. We take the 3 atoms to have the same mass m. 
The 3-body recombination rate R{pi,p2,P3) is a function of the momenta of the three 
incoming atoms. Galilean invariance implies that R does not depend on the total momentum 
Ptot — Pi + P2 + P3- It can therefore be expressed as a function -R(pi2, ^3,12) of a pair of 
Jacobi momenta pi2 = Pi — P2 and ^3,12 = Ps — ^{Pi +^2)- It is convenient to parameterize 
the Jacobi momenta in terms of 6 orthogonal variables: the collision energy E = (3p^2 + 
^2)/(12m), four Jacobi angles giving the orientations of the unit vectors pi2 and ^3,12, 
and a hyperanglc = arctan(\/3pi2/(2p3.i2))- The integration clement in these variables is 

d^Pi d^P2 d^Ps = VSm^E^dE sm'^{2a3)da3 dVLi2 dVt^^i2 d^ptot- (2) 

We will refer to the average of a quantity over the Jacobi angles and over the hyperangle as 
the hyperangular average. We denote the hyperangular average of the 3-body recombination 
rate i?(pi2, ^3,12) by K{E). 

In a gas of atoms with number density ua, the rate of decrease in the number density 
due to 3-body recombination defines an experimentally measurable loss rate constant L3: 

^riA = -L3 n\ . (3) 

The event rate constant a for 3-body recombination is defined so that the number of recom- 
bination events per unit volume and per unit time is an\. If the number of atoms lost from 
the system per recombination event is niost, the rate constant is L3 — niost ck- The binding 
energy of the dimer is released through the kinetic energies of the recoiling atom and dimer. 
If their kinetic energies are large enough that the atom and dimer both escape from the 
system and if the number density is low enough that subsequent collisions do not prevent 
them from escaping from the system, then niost — 3. On the other hand, if they both remain 
in the system and if we regard the dimer as a distinct chemical species with its own number 
density no, then niost = 2. 



Definitive zero-range predictions at temperatures small compared to the dimer binding energy are available 
for resonant dimer relaxation, which is an important loss process for atom-dimer mixtures [25] . 
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In an ensemble of identical bosons with number density ua in thermal equilibrium at 
temperature T, the event rate per unit volume and per unit time is an% where 



^ / R{Pi2, P3,i2) n(Ei)n(E2)n(E3) d^pid^p2d^P3 . . 

^ 3\ fn{E,)n{E2)n{Es)d^Pid^P2d^Ps 

and n{E) is the Bose- Einstein distribution: n{E) = [e^'^"'^)/^'^^'^) — The chemical 

potential fj, is determined by the condition J n{Ei)d^pi/{2-K)^ = ua- The factor of 1/3! 
in Eq. (4) accounts for the indistinguishable nature of the three identical particles. The 
denominator in Eq. (4) can be written 6 [J n{Ei) rf'^pi]'^ = 6(27r)^n^. 

The critical temperature Tc for Bose-Einstcin condensation is given by ksTc ~ 
S.Sh^TiA^ /m. If T is significantly larger than the product of the three Bose- Einstein 
distributions in Eq. (4) can be approximated by a Boltzmann distribution: 

n{Ei)n{E2)n{Es) ^ e^'^/^'^-^) exp(-^tot/(A;BT)), (5) 

where E'tot = Ei + E2 + E3 is the total energy of the three atoms. The total energy can 
be expressed as the sum of the center-of-mass energy and the collision energy E: E^ot = 
Ptot/(6^) + E. The Gaussian integral over ptot cancels between the numerator and the 
denominator of Eq. (4). The effect of the integrals over the hyperangle and the Jacobi 
angles is to replace i?(pi2,P3,i2) in the numerator by its hyperangular average K{E). Thus 
if T is significantly larger than the event rate constant a{T) in Eq. (4) reduces to a 
Boltzmann average over the collision energy E: 

J^dEE^e-^/^'^'^^^KiE) 
""^^^ 6 J^^ dE E^ e-^/C^BT) ■ (6) 

The integral in the denominator can be evaluated analytically to give 2{kBT)^. We will 
refer to temperatures and densities for which this approximation is valid as the Boltzmann 
region. Note that the Boltzmann region can extend to arbitrarily low temperatures if the 
number density ua is arbitrarily small. At T = 0, the approximation in Eq. (6) reduces 
to a(0) = K{E = 0)/6. However this applies only if the temperature is small enough that 
K{E) can be approximated by K{0) but still in the Boltzmann region. At T = 0, the atoms 
are in a Bose-Einstein condensate. The correct result for the rate constant at T = is 
therefore a{0) K{E — 0)/36. The extra factor of 1/3! comes from the 3 identical bosons 
being in the same quantum state. 

If the scattering length a is positive and large compared to the range of the interaction, 
one of the dimers that can be produced by the recombination process is the shallow dimer. 
Its binding energy in the zero-range limit is 

Ed = hy{ma^) . (7) 

If there are deep dimers, they can be produced by the recombination process for either sign 
of a. The recombination rate can be decomposed into the contribution from the shallow 
dimer and the sum of the contributions from all the deep dimers: 

K{E) = Xshallow(£^) + K^eeAE) ■ (8) 
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The 3-body recombination rate into the shallow dimer can be further decomposed into 
contributions from the channels in which the total orbital angular momentum of the three 
atoms has definite quantum number J: 

oo 

i^shaIIow(£^)= J^i^^-'H^)- (9) 

J=0 

The threshold behavior for each of the angular momentum contributions to i^^ghaiiow(£') 
follows from a generalization of Wigner 's threshold law [26]: K^-^\E) ~ E^' , where Aq = 0, 
Ai = 3, and Aj = J for J > 2. At the scattering threshold E = only the J = term is 
nonzero. 

The hyperangular average -f^shaiiow(-E') of the 3-body recombination rate is related in a 
simple way to the dimer-breakup cross section for the scattering of an atom and the shallow 
dimer: 

lQ2^7:n\En + E) 

-f^ shallow I -C'j = ^2^2 CTbreakupl-C/J- \^^) 

Thus in the Boltzmann region, the contribution to the event rate constant in Eq. (6) from 
3-body recombination into the shallow dimer is determined by the dimer-breakup cross 
section o"breakup(-E)- Usiug the Optical Theorem, o-brcakup(-£') can be determined from the 
phase shifts for elastic atom-dimer scattering. We will show in Section VI that in the case 
of a large scattering length, K^^epiE) is determined by the same universal scaling functions 
that determine -f^shaiiow(-E')- Thus in the Boltzmann region, the rate constant for 3-body 
recombination is completely determined by the elastic atom-dimer scattering amplitude. 



III. UNIVERSAL SCALING BEHAVIOR 

In this section, we consider atoms with a large positive scattering length and no deep 
dimers, so the only diatomic molecule is the shallow dimer whose binding energy is given by 
Eq. (7). The total 3-body recombination rate K{E) is therefore equal to -ft'shaiiow(-E'), which 
is related to the dimer-breakup cross section by Eq. (10). We summarize the constraints of 
universality on the 3-body recombination rate for this case. 



A. Angular momentum and hyperangular decompositions 

The zero-range results for scattering rates, such as -f^shaiiow(-E')) ^-re functions of the col- 
lision energy E^ the scattering length a, and the Efimov parameter defined by Eq. (1). 
An alternative Efimov parameter that is particularly convenient when considering 3-body 
recombination is a*o, which is one of the values of the scattering length at which Kshaaov,{E) 
vanishes at £■ = 0. This Efimov parameter differs from by a multiphcative factor that 
is known only to two digits of accuracy [2] : 

a*o~0.32«;;\ (11) 

The dependence of scattering rates on the Efimov parameter is strongly constrained by the 
discrete scale invariance. Their dependence on the energy E can be conveniently expressed 
in terms of universal functions of a dimensionless scaling variable x defined by 

x^{ma''E/h^Y'\ (12) 
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Atom-dimer scattering states in the center-of-mass frame can be labelled by the relative 
momentum k between the atom and the dimer or, equivalently, by the energy E relative 
to the 3-atom threshold and the angular momentum quantum numbers J and M. The S 
matrix elements for atom-dimer elastic scattering are functions of E that are independent 
of M and diagonal in J. They can be expressed in terms of the phase shifts 5^fl,{E) for 
atom-dimer elastic scattering: 

SfnME) = e"^-^^^ (13) 

The phase shifts 5^^1{E) are real below the 3-atom threshold and have a positive imaginary 

part for £' > 0. 

Three-atom scattering states in the center-of-mass frame can be labelled by the Jacobi 
momenta pi2 and P3,i2- One can define subsystem angular momenta associated with ro- 
tations of the unit vectors pi2 and P3,i2- We denote the associated angular momentum 
quantum numbers by i^, and by iy, rriy, respectively. An equivalent set of quantum 
numbers are £x, and the total orbital angular momentum quantum numbers J and M. 
The possible values of J are integers ranging from \ix — £y\ to + iy. The range of the 
hyper angle is from to ^tt. Since this range is compact, we can expand functions of 
in terms of orthogonal functions relative to the weight factor sin^(2Q;3) that are labelled by 
an index = 0, 1,2, . . .. A convenient set of variables for the 3-atom scattering states is 
the collision energy E and five discrete variables: J, M, i^, iy, and the hyperangular index 
713 . Since ix, iy, and na are a denumerable set, we will denote them collectively by a single 

integer n that takes values 3, 4,5, The reason for choosing this pecuhar set of values 

will become clear in Section III C. The S-matrix elements for 3-body recombination can be 
expressed as functions S^^^^]^ y^j^{E) of the collision energy with index n. The contribution 
from angular momentum J to the hyperangular average of the 3-body recombination rate 
can be expressed as 

n=3 

The hyperangular average is implemented by the sum over n. The unitarity of the S-matrix 
in the angular momentum J sector implies 



n=3 



This can be used together with Eq. (13) to express the 3-body recombination rate in Eq. (14) 
in terms of the phase shifts (^^(-E") for elastic atom-dimer scattering: 

m.'^E'^ \ J 



B. Scaling behavior for J > 1 

In the sector with angular momentum quantum number J > 1, the scattering rates do 
not depend on k*. The constraints of universality are therefore particularly simple for J > 1. 



7 



The atom-dimer phase shifts (^^^(i?) defined by Eq. (13) are universal functions of the scahng 
variable x defined in Eq. (12). Below the 3-atom threshold, these scaling functions are real. 
Above the 3-atom threshold, they are complex. 

The contribution to the 3-body recombination rate into the shallow dimer from angular 
momentum J has the form 

^(,) ^ lUV37r\2J+l)fj{x) ha^ 

where fj{x) is a real- valued scaling function: 

fj{x) = 1 - exp ( -Almd^iliE)) . (18) 

As X 0, the leading powers of x are determined by Wigner's threshold law [26]: fj{x) ~ 
^2Aj+4^ where Ai = 3 and Aj = J for J > 2. 



C. Discrete scaling behavior for J = 

The zero-range results for J = are more intricate than for J > 1, because they depend 
not only on E and a, but also on or, equivalently, a^o- General scaling formulas for scat- 
tering rates in the 3-atom sector can be derived from Efimov's Radial Law. The Radial Law 
expresses the S-matrix elements for low-energy scattering processes in the J = channel of 
the 3-atom sector in terms of universal functions of the scaling variable x defined in Eq. (12). 
The S-matrix elements for atom-dimer elastic scattering and for 3-body recombination from 
3-atom states labelled (J = 0, n) have the form [2]: 

where the angle 9^q is 

= soln(a/a*o). (20) 

The functions Sij{x) are entries of an infinite-dimensional symmetric unitary matrix that 
depends only on the scaling variable x defined in Eq. (12). The unitarity of the infinite- 
dimensional matrix s can be expressed as 

oo 

^li^lj + *2i*2i + ^ S*^iSrij = ^ij- (21) 
n=3 

This condition implies the unitarity of the physical S-matrix for the J = sector. For 
example, one can use Eq. (21) to verify that the unitarity condition in Eq. (15) for J = is 
automatically satisfied if the S-matrix elements are given by the expressions in Eqs. (19). 

Efimov's Radial Law has a simple interpretation in terms of the adiabatic hyperspherical 
representation of the 3-atom Schrodinger equation. The hyperradius R — [(rj2 + '"23 + 
r|i)/3]^/2 jg the root -mean-square separation of the three atoms. We take the scattering 
length a to be much larger than the absolute value of the effective range r^: a ^ \rs\. In 
this case, there are four important regions of R for 3-atom configurations with energies near 
the scattering threshold: 
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• the asymptotic region a, 

• the scaling region R a, 

• the scale-invariant region \rs\ <^ R ^ a, 

• the short- distance region R 

Efimov's radial law reflects the fact that there is only one adiabatic hyperspherical potential 
that is attractive in the scale-invariant region of R. It is only through this adiabatic hyper- 
spherical potential that 3-atom configurations with energies near the scattering threshold 
can reach the short- distance region of R. The outgoing and incoming asymptotic states with 
respect to the scaling region are 

1. incoming and outgoing hyperradial waves in the scale- invariant region \rs\ <^ i? ^ a 
of the lowest adiabatic hyperspherical potential. We label these asymptotic states by 
the index 1. 

2. outgoing or incoming atom-dimer scattering states with angular momentum quantum 
number J = 0. We label these asymptotic states AD or simply by the index 2. 

3. outgoing or incoming 3-atom scattering states with J = 0. We label these asymptotic 
states AAA by the index n = 3, 4, 5, . . . that specifies i^, ^y, and Us. 

The evolution of the 3-atom wavcfunction through the scaling region is described by a 
unitary matrix Sij whose indices correspond to the asymptotic states relative to the scaling 
region that were enumerated above. Time-reversal invariance implies that the matrix s is 
symmetric. 

The expressions for the S-matrix elements in Eq. (19) have simple interpretations. The 
factor e^*^*" in Eqs. (19) is the phase shift due to reflection of a hyperradial wave from the 
short- distance region. The S-matrix elements in Eqs. (19) can be expanded as a power 
series in e^*^*° with each term corresponding to a different pathway between the incoming 
and outgoing scattering states. The term with the factor (e^*^*o)" is the amplitude for a 
pathway that includes n reflections from the short-distance region. In the S-matrix element 
^AD^ADi-^) ^'i- (19a), the leading term S22 is the amplitude for the incoming atom-dimer 
scattering state to be reflected from the scaling region. The second term S2ie^*^*"si2 is 
the amplitude for the atom-dimer scattering state to be transmitted through the scaling 
region to an incoming hyperradial wave, which is reflected from the short-distance region 
into an outgoing hyperradial wave, which is then transmitted through the scaling region 
to an outgoing atom-dimer scattering state. The factors of e^*^*°sii from expanding the 
denominator are the amplitudes for the hyperradial wave to be reflected from the short- 
distance region and then reflected from the scaling region. 

The contribution to the hypcrangular average of the 3-body recombination rate from 
angular momentum J = can be obtained by using the unitarity condition in Eq. (15) to 
eliminate the sum over n in Eq. (14) and then inserting the expression for the S-matrix 
element in Eq. (19a): 



suixfe'^'*'^ 



1 - sii(x)e2^^*o 




(22) 
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D. Analytic results near threshold 



Macek, Ovchinnikov, and Gasaneo have derived some remarkable analytic results for the 
3-body problem of identical bosons with large scattering length and energies near the 3- 
atom threshold [27, 28]. In Ref. [27], they derived an analytic result for the S-wave phase 
shift ^^'^'^''(0) for atom-dimer elastic scattering at the breakup threshold E = 0. They 
obtained analytic results for the values Sii(O) and Si2(0) of the universal scaling functions at 
this threshold. In Ref. [28], they extended their analytic solution for S^^^^\e) to the first 
nontrivial order in E. They used this result to derive an analytic expression for the 3-body 
recombination rate i^^shaiiow(O) at the 3-atom threshold. 

An analytic expression for the atom-dimer S-wave phase shift 6^^jj{E) near the dimer- 
breakup threshold is given in Ref. [28]. Using the first equality in Eq. (54) of Ref. [28] 
(to avoid typographical errors in the last line of Eq. (54)), the S-matrix element at energy 
E — x^Ed can be expressed in the form 

exp {2iS%{E)) = ^ ^'^";!te-2^ (.) + e-^-°e^'^(^) + e--A(x) (e^^^(^) - l)] , (23) 
where 

A{x) = 6o{x)-soHa/Ro). (24) 

The variable Rq is a matching point at short distances where the 3-body wavefunction 
vanishes. (The function A(x) was denoted by A(i?o) in Ref. [28].) The function A{x) and 

the energy dependent phases ^00(2^) and Sq{x) are universal functions of the scaling variable 
X. By replacing the expression inside the absolute value signs in Eq. (22) by the phase shift 
in Eq. (23), the 3-body recombination rate can be expressed in the form 

K^^\E) = 288V37r^^ ( 1 - (^^0) ^^^^ ^(") ^ . (25) 

V 2sinh(7rso)y sinh^(7rso) + cos2 A(a;) m 

To determine the 3-body recombination rate at threshold, we need the limiting behavior 
of the functions A{x) and So{x) as x — > 0. The values of these phases at a; = were 
calculated numerically in Ref. [27]: 

S^(0) = 1.736, (26a) 
5o(0) = 1.588. (26b) 

The leading term in ^4 (a;) as a; ^ 0, which is proportional to x'^, was calculated analytically 
in Ref. 



Ai,) ^ ^(ll^ML,., (27) 

3v3 sinh(7rso) 

The recombination rate in Eq. (25) vanishes at the threshold if A(0) = 0. Thus the expres- 
sion for A(x) in Eq. (24) can be written 

A{x) = 5o{x) - 5o{0) - So ln(a/a,o) , (28) 

where a*o is a value of the scattering length for which the 3-body recombination rate vanishes 
at threshold. The expression for the S-wave atom-dimer phase shift in Eq. (23) reduces at 
the 3-atom threshold to 

S%{E = 0) = 5oo(0) + arctan (tanh(7rso) tan^^o) ■ (29) 
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The expression given in Ref. [27] is equivalent but less compact. The analytic result of Macek, 
Ovchinnikov, and Gasaneo for the 3-body recombination rate at the 3-atom threshold [28] 
can be obtained by taking the limit a; — > in Eq. (25): 

^(0) _ _ 7687r^(47r - 3v^) sin^[go ln(a/a,o)] ha^ ^^^^ 
sinh^(7rso) + cos^[soln(a/a*o)] 

This result was also derived independently by Petrov [29]. For fixed a, the rate in Eq. (30) 
is a log-periodic function of a*o that oscillates between zero and a maximum value 

^max = QC^^ha'^/ m, (31) 

where C^ax is 



_ 1287r^(47r-3^) 
smh (ttso) 

Its numerical value is Cmax ~ 67.1. The expression in Eq. (30) has zeroes when a is 
(e'^/'*")"' a*o, where n is an integer. The maxima of K^'^\0) in Eq. (33) occur when a is 
(e'^/^o)" 14.3 0*0- Since sinh^(7rso) !^ 139 is so large, the expression for K^^\0) in Eq. (30) 
can be approximated with an error of less than 1% of QC^sJ^a^ / m by 

K^^\E = 0) fti 6C„,axSin2[soln(a/a,o)] fia^/m. (33) 

This approximate functional form of the rate constant was deduced independently in 
Refs. [6-8]. 

By comparing the expression for the S-matrix element in Eq. (23) with the Efimov's 
Radial Law in Eq. (19a), we can determine the universal scaling functions sii(x), 512(2;), 
and 522(2;) in the approximation of Ref. [28]: 



si2{x) ^ e^f'^°°('^)-^°('')+^°(°^Vl - 



A{x) 



.1/2 



(34b) 



2 sinh(7rso)_ 

S22{x) ^ e^^ioo{x)^-^2^so [1 + e"^«A(a;)] . (34c) 

The limiting behavior of the functions 5oo{x)., 5q{x)^ and A{x) as a; — > 0, is given in 
Eqs. (26) and (27). At a; = 0, we obtain sii(O) = -e-^^^^" and Si2(0) = e'^^^^\l - e-^^'" f'^ . 
These differ from the values deduced in Ref. [27] because a phase e^*''"''^-* has been absorbed 
into e^*^*°. They also differ by a minus sign in Sii(O). For small x, the leading terms in 
512(2;) — 512(0) and 522(2;) — ^22(0) scale like x^. The leading terms in ^11(2;) — Sii(O) scale 
like a higher power of x. 



IV. UNIVERSAL SCALING FUNCTIONS 

In this Section, we calculate the universal scaling functions associated with atom-dimer 
elastic scattering. These functions are then used to calculate the hyperangular average of 
the 3-body recombination rate as a function of the collision energy. 
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A. STM Equation 



Three-body observables for systems with a large scattering length can be calculated by 
solving an integral equation called the Skorniakov-Ter-Martirosian (STM) equation [30]. In 
momentum space, the STM equation can be expressed in the form [2] 



mE — {p^ + p ■ k + k'^) + ie 

d^q Stt A{q, k; E) 



I 



{2ny mE - {p"^ + p ■ q + q"^) + ie -i / a + ,/^^^^^eTWJ^^ 



.(35) 



The solution is, up to a normalization constant, the amplitude for an atom with momentum 
p and energy p^/ (2m) and a pair of atoms with total momentum —p and total energy 
E — / (2m) to evolve into an atom with momentum k and energy A;^/ (2m) and a pair of 
atoms with total momentum —k and total energy E — k'^/{2m). The STM amplitude can 
be resolved into the contributions from channels with orbital angular momentum quantum 
number J: ^ 

Aj{p, k;E) = lJ ^dy Pj{y)A{p, k; E) , (36) 

where y — p • k/pk and Pj{y) is a Legendre polynomial. 

In a specific angular momentum channel, the STM equation in Eq. (35) reduces to an 
integral equation with a single integration variable q. In the case J = 0, the behavior of the 
integral is sufficiently singular at large q that it is necessary to impose an ultraviolet cutoff 
g < A, where A is much greater than p, k and {m\E\Y/'^. The STM equation for J = 
reduces to 

. , , „ . X Stt , + pfc + /c^ - mE — ie 

TT Jo P P^ -pq + q^ -mE -i€ -l/a + y/Sq^/A - mE - ie 

It can be shown that changing the ultraviolet cutoff A corresponds to changing the Efimov 
parameter [31-33]. The values of A and differ simply by a multiplicative constant: 

ln(A/K^) ^ 1.74 mod (tt/sq) . (38) 

The solutions of Eq. (37) are log-periodic functions of A. If A is increased by a factor of 
g7r/so ~ 22.7, it corresponds to the same value of k^. The S-wave atom-dimer phase shift 
6^^jj{E) can be obtained from the amplitude with both momentum arguments p and k set 
equal to ks = {lm{ED + E)f/'^: 

37r 

Ao{kE, kE-, E, A) = .0),^, ., ■ (39) 

kE cot d\j-){E) — ikE 

This phase shift depends not only on the parameter a but also on the parameter through 
the dependence of the amplitude on A. 
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In the case J > 1, no ultraviolet cutoff is required and the partial wave projection of the 
STM equation leads to (see, e.g. Ref. [34]) 



Trio V pq / -1/a + ^3^74 - mE - ie 



where 

Qj{z) dx^ (41) 

2J_i z-x 

is a Legendre function of the second kind. The atom-dimer phase shift 5^ad{E) can be 
obtained from 

37r 

MkE^E-.E)^- , for J>1. (42) 

kE cot o\l){E) — ikE 

This phase shift is a function of the dimensionless variable ma'^E/h^. 



B. Universal scaling functions for J > 1 

We calculate the atom-dimer phase shifts S^^^{E) for J — 1,2, ... ,6 as functions of the 
energy E from lO^'^Ejy to IO^Ed. The scaling variable x defined in Eq. (12) then ranges 
from 0.1 to 10. For each energy E, we solve the STM equation in Eq. (40) and determine the 
phase shift 6^/^{E) using Eq. (39). The universal scaling function fj{x) defined in Eq. (18) 
is then determined from the imaginary part of the phase shift. Our numerical method loses 
accuracy for x smaller than values that range from 0.1 for J = 2 to x = 0.45 for J — 6. To 
determine fj{x) for smaller values of x, we fit our results for the lowest 10 accurate points 
to the form 

fj(x) = ajx^^-'+^ + hjx^^^^\ (43) 

where Ai = 3 and Aj = J for J > 2. We use the formula in Eq. (43) to extrapolate to small 
values of x where our numerical method loses accuracy. Our results are shown in Fig. 1. For 
J >2, fj{x) is a decreasing function of J. This pattern is broken by fi{x), which is smaller 
than f2{x) and fsix) at small x but eventually becomes larger than f2{x) around x — 9. 



C. Universal scaling functions for J — 

We calculate the atom-dimer phase shift 6^^1){E) by solving the STM equation for J = 
with an ultraviolet cutoff A, which is given in Eq. (37), for energies E ranging from 10~^Eo 
to IO'^Ee). We choose values of the ultraviolet cutoff A that are equally spaced on a log 
scale and cover three quarters of a discrete scaling cycle from Aq to 10.4Ao. The smallest 
cutoff Ao is chosen to be much greater than 1/a and (■m\E\/h'^y^'^ and large enough that 
the calculated phaseshift is unchanged if A is increased by a factor of 22.7. For each energy 
E, we calculate the phase shift as a function of log(aA) and fit it to the universal formula 

/ .^(j=o).px\ ... Si2(a;)2exp[2isoln(a/a*o)] .... 

exp 2id\r, (E)] = S22ix) + -—^ f— — t——, -T (44) 

V ^ V ^ ^ 1 - sii X exp 2zso In a/a*o 
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FIG. 1: The universal scaling functions fj{x) = 1 — exp{—4Im.6 j^^{E)) for J = 1, . . . , 6 as functions 
of the scaling variable x. 

to determine the universal scaling functions Si2(x), S22{x), and Sii(x). The sign ambiguity 
in Si2(x) is resolved by using continuity together with the analytic result for Si2(0) given by 
Eq. (34b). 

The numerical results for the modulus and the phase a.TgSij{x) of each of these 

functions are shown in Fig. 2. The modulus |si2(2^)| remains close to 1 for x < 1 but it 
decreases to 0.59 at x = 10. The moduli |sii(x)| and |s22(a;)| both have the tiny value 
0.0018 at a; = and they remain tiny for x < 1. While |s22(a^)| increases to 0.37 at a; = 10, 
increases much more slowly to 0.035. 

D. Three-body recombination rates 

In this subsection, we use the universal scaling functions calculated in the previous sec- 
tions to calculate the 3-body recombination rate into the shallow dimer. We show the 
contributions from angular momentum J = through 6 as functions of the collision energy. 

In Eq. (22), the J = contribution K^^\E) to the 3-body recombination rate is expressed 
in terms of the universal scaling functions sii(x), Si2{x), and ^22(2;). In Fig. 3, it is shown 
as a function of the collision energy E for 6 values of O^^q: 0, ^vr, ^tt, |7r, ^tt, and ||7r. For 
Q*o = ^TT, K^^\E) at E = has the maximum possible value i^max given in Eq. (31) and it 
decreases monotonically with E. For = or vr, K^^\E) vanishes at E = 0, increases to 
a maximum value of 0.0143 -ft'max aX E = 4.61 E'o, and then decreases. For between |7r 
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FIG. 2: The universal scaling functions si2{x) (upper panels), S22{x) (middle panels), and sii(x) 
(lower panels) as functions of x. The modulus |sij(x)| and the phase argSij(x) of the function are 
shown in the left panel and the right panel, respectively. 
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FIG. 3: (Color online) The J = contribution to the 3-body recombination rate into the shallow 
dimer for the case in which there are no deep dimers. The hyperangular average K^'^\E) is shown 
as a function of the collision energy E for 6 values of O^q: ^tt (upper solid line), ^tt and j^tt (upper 
and lower dashed lines), |^7r and ^vr (upper and lower dash-dotted lines), or vr (lower solid line). 

and vr, K^^\E) has its maximum at a nonzero energy. As 9^:0 increases from ^vr to vr, the 
local minimum at £^ = becomes increasingly deep and the maximum moves outward from 
to 4.61 Ed. For 9,o between and 0.124??, K^^\E) has a local minimum at a nonzero 
energy E. As increases from to 0.1247r, the local minimum moves outward from to 
about 8 Ed and its depth decreases until it becomes just an inflection point. For 9^,^ between 
0.124-7? and |7r, K^^^E) decreases monotonically with E. 

At small energies, the dimensionless recombination rate K'^^\E)/K^^^ varies by many 
orders of magnitude as a=KO ranges over a complete discrete scaling cycle. Even at ii^ = Ed, 
it varies by more than 3 orders of magnitude. The variations with a^^Q become smaller at 
larger energy. If there were no relevant length scales at high energy, dimensional analysis 
would imply that K'^^\E) should be proportional to E'"^ at large E. Our results at the 
largest values of E are compatible with the approach to this simple scaling behavior. The 
differences between K^^\E) for different values of o^kq seem to decrease as a higher power 
of E. Our results at the largest values of E are compatible with their approach to the 
asymptotic behavior E~^^'^. 

In Eq. (17), the contribution K^-^^E) to the 3-body recombination rate from angular 
momentum J is expressed in terms of the universal scaling functions fj{x). In Fig. 4, the 
rates K^'^\E) for J from 1 to 6 are shown as functions of the collision energy E. Over the 
range from to 100 Ed, the largest of these contributions is the J = 2 term. It reaches its 
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FIG. 4: (Color online) The contributions to the 3-body recombination rate into the shallow dimer 
from angular momenta J = 1, ... ,6. The hyperangular averages K^'^\E) are shown as a function 
of the collision energy E. 

maximum value of 0.214 -ft'max aX E = 4.18 Ed. The total contribution from J = 1, 3, 4, 5, 
and 6 is less than 10% of the J = 2 contribution ii E < 0.78 Eo- Thus the sum of the terms 
for J = and 2 is a good approximation to the recombination rate if £^ < 0.78 Ejj. The 
sum of the contributions from J = 5 and 6 is less than 10% of the sum of the contribution 
from J = 1 through 4 if E' < 32 E^. Thus the sum of the terms for J = through 4 is 
a good approximation to the recombination rate if < 32 Ed. Truncations of the partial 
wave expansion quickly become inaccurate at higher energies. 

V. APPLICATION TO ^HE ATOMS 

The results from the previous section can be applied directly to systems of ^He atoms. ^ 
There are several modern potentials that are believed to describe the interactions of ^He 
atoms with low energy accurately. They all support a single weakly-bound diatomic molecule 
(or dimer). The effective ranges for all these potentials are approximately 14 oq, 

which is a little larger than the van der Waals length scale {mCfs/h^Y'^ = 10.2 gq. The 
scattering lengths are larger by about a factor of ten, and they vary among the potentials. 
The large scattering length, a 3> r^, implies that ^He atoms have universal properties 



^ A convenient conversion constant for "^He atoms is h^/m = 43.2788 K Qq 
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that are determined by a and the Efimov parameter k^,. First-order range corrections to 
the universal results are estimated to be suppressed by roughly r^/a ii E ^ En and by 
Tg {mE /Tl^Y^"^ = xTs/a for E ^ En- This makes ^He atoms an ideal system to illustrate 
the universal properties of atoms with large scattering length. A thorough analysis of the 
universal properties of ^He atoms has been presented by Braaten and Hammer [35] . Various 
3-body observables for ^He atoms have been calculated by Platter and Phillips to next-to- 
next-to-leading order in the expansion in powers of r^/a [13]. The agreement with numerical 
results from exact solutions of the 3-body Schrodinger equation [36, 37] is impressive. 

The 3-body recombination rate for ^He atoms was first calculated by Suno, Esry, Greene, 
and Burke [38]. They solved the 3-body Schrodinger equation for ^He atoms interacting 
through the HFD-B3-FCI1 potential [39]. They calculated /rshaiiow(-£') as a function of the 
collision energy from the threshold to 10 mK [38], including the contributions with angular 
momentum quantum number J = 0, 1, 2, and 3. Their results are shown as dots in Fig. 5. 
Shepard has calculated the 3-body recombination rates of ^He atoms for the HFD-B3-FCI1 
and other ^He potentials using separable potentials that reproduce the corresponding two- 
body phase shifts together with a 3-body force that is adjusted to fit the binding energy of 
the excited trimer [23]. 
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FIG. 5: (Color online) The 3-body recombination rates K^'^\E) for ^He atoms (in units of cm^/s) 
as functions of the collision energy E (in units of mK) for J = through 6. The solid lines are 
our zero-range results for a = 164.5 oq and a=KO = 143.1 oq. The dots are the results obtained in 
Ref. [38] by solving the 3-body Schrodinger equation for the HFD-B3-FCI1 potential. 

To apply the zero-range predictions for the 3-body recombination rate to ^He atoms 
interacting through the HFD-B3-FCI1 potential, we must determine the scattering length 
a and the 3-body parameter a*o for that potential. The scattering length for the HFD-B3- 
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FCIl potential is a = 172 oq and the binding energy of the shallow dimer is En — 1.600 
mK [40]. The zero-range predictions for low-energy 3-body observables are most accurate 
if the dimer bound state pole is at the correct position [35]. Therefore, we use the dimer 
binding energy as the 2-body input instead of a. Using the zero-range expression for the 
dimer binding energy in Eq. (7), we obtain the scattering length 

a"" = 164.5 ao . (45) 

We can determine the value of a*o for ^He atoms interacting through the HFD-B3-FCI1 
potential from the value of the 3-body recombination loss rate at threshold: -ft'shaiiow(O) = 
7.10 X 10'^^ cm^/s. Inserting a^'^ = 164.5 Oq in Eq. (30) and solving for a^Q, we obtain 

= 143.1 ao. (46) 

The result for af^ is approximately equal to 0.870 a^^. The near equality between af^ 
and a^^ reflects the fact that -ft'shaiiow(O) for ^He is much smaller than the maximum value 
6Cmax^('^^'^)^/'^ = 3.67 X 10~^^ cm^/s allowed by universality, which implies that a/a*o is 
close to a zero of the sine function in the numerator of Eq. (30). The value of a*o for ^He can 
also be determined from the binding energy of the excited "^He trimer, which is E^^ = 2.62 
mK for the HFD-B3-FCI1 potential [40]. The resulting value, a*o ~ 146 oq, is consistent 
with the value in Eq. (46). Its accuracy is limited by the accuracy to which the numerical 
value of the constant in Eq. (11) is known. 

The zero-range predictions for K^'^\E) are given in Eq. (22) for J = and in Eq. (17) 
for J > 1. The only information about the HFD-B3-FCI1 potential that we use are the 
values of a and a*o given in Eqs. (45) and (46). Our results for K^''\E), J = 0, 1, . . . , 6, 
are shown as solid lines in Fig. 5. The J = term dominates at low energies, but has a 
dramatic dip near 2 mK, where it is smaller than its value at ii^ = by about a factor of 30. 
The J = 2 term takes over as the dominant term when the energy exceeds about 0.2 mK. 
When it reaches its maximum value at an energy near 4 mK, other partial waves begin to 
give important contributions. 

Our results can be compared with those of Ref. [38] for J = 0, 1, 2, and 3, which are shown 
as dots in Fig. 5. The qualitative agreement is excellent. The deep local minimum for J = 
and the peculiar shape for J = 1 are both reproduced by the zero-range results. There is 
also good quantitative agreement for J = 0, 1, and 2. At = Ed — 1.6 mK, the fractional 
discrepancies for J = 0, 1, and 2 are 18.5%, 14.4%, and 20.3%, respectively. This is a 
little larger than the expected fractional error from range corrections, which is r^/a = 8.5%. 
At £" = 10 mK, which is the highest energy for which the 3-body recombination rate was 
calculated in Ref. [38], the fractional discrepancies for J = 0, 1, and 2 are 16.5%, 6.0%, and 
18.6%. They are consistent with the expected fractional error from range corrections, which 
is XTs/a = 21%. For J = 3, the agreement between our results and those of Ref. [38] is not as 
good. The fractional discrepancy varies from +77% at 0.1 mK to —36% at 10 mK. Shepard 
has also calculated the contributions to the 3-body recombination rate for the HFD-B3-FCI1 
potential from partial waves up to J = 3 [23] . His results for J = 3 are similar to ours and 
show the same discrepancy with Ref. [38]. The agreement between our results for J = 3 and 
those of Ref. [23] gives us confidence in the accuracy of our numerical calculations. We do 
not understand the origin of the discrepancy with the results of Ref. [38]. 

The inclusion of higher order corrections in r^/a should improve the agreement with the 
results of Ref. [38] at all energies. If the first order corrections in Vg/ a are included, the errors 
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should decrease to roughly (vs/a)^ = 0.7% for E ^ Er> = 1.6 mK and to roughly {xrs/a)^ at 
higher energies. A first calculation of the range corrections to the threshold recombination 
rate i^shaiiow(O) was carried out in Ref. [41]. 



VI. EFFECTS OF DEEP DIMERS 

In this section and in the subsequent section, we consider atoms that have deep dimers. 

In the scaling (or zero-range) limit, the cumulative effect of all the deep dimers on Efimov 
physics can be taken into account through one additional parameter rj^ [15]. The atom-dimer 
phase shift and the hyperangular average of the 3-body recombination rate are determined 
by the same universal scaling functions as for 77* = 0. 



A. Generalization of Efimov's Radial Law 

The existence of deep dimers implies that there are states in the 3-atom sector with 
energies near the 3-atom threshold that consist of an atom and a deep dimer with large 
kinetic energies. We will denote these states by the symbol AD. These states provide 
inelastic atom-dimer scattering channels and additional 3-body recombination channels. In 
the scaling limit, the only contributions to the recombination rate into deep dimers are from 
the J = channels. The contributions from J > 1 are suppressed by powers of r^/a. The 
hyperangular average of the inclusive 3-body recombination rate into deep dimers can be 
expressed as 

l il\/3Ti^h^ ^ ^ I „(J=0,n) 2 



n=3 An 

where 'S'^aa Ar)(-^) the S-matrix element for the transition from a 3-atom scattering state 
labelled by (J — 0,n) to a specific state AD consisting of an atom and deep dimer. The 
hyperangular average is implemented in Eq. (47) by the sum over n. 

The states AD also have an effect on the amplitudes for atom-dimer elastic scattering and 
for 3-body recombination into the shallow dimer. These effects are strongly constrained by 
the fact that the states AD can be accessed only through a single adiabatic hyperspherical 
channel: the J — channel whose potential is attractive in the scale-invariant region of 
the hyperradius R. If there are deep dimers, an incoming hyperradial wave need not be 
totally reflected from the short-distance region of R, because it can ultimately emerge as a 
scattering state consisting of an atom and a deep dimer. The reflected hyperradial wave will 
therefore differ from the incident hyperradial wave not only by a phase shift e^'^*° but also 
by a decrease in its amplitude by a factor e"^''*. The parameter r)^, which was introduced 
in Ref. [14], determines the widths of Eflmov trimers as well as all other effects of the deep 
dimers on low-energy 3-atom scattering processes. If the zero-range result for an amplitude 
for the case of no deep dimers is known as an analytic function of a*o, the corresponding 
amplitude for a system with deep dimers can be obtained without any additional calculation 
simply by making the substitution 

In a*o — >■ In a^o - iv*/so ■ (48) 
The generalization of Efimov's Radial Laws in Eqs. (19) to the case of atoms with deep 
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dimers are 



c(^=0,n) , . . 52l(x)5ln(x)e^^^-°-^^- 

Thus these S- matrix elements are determined by the same universal scaling functions Sjj (x) 
as in the case ry* = with no deep dimers. 

The S-matrix elements for transitions from an atom-dimcr scattering state or from a 
3-atom scattering state into a specific state AD consisting of an atom and a deep dimer 
depend on details of physics at short distances. However as pointed out in Refs. [14, 15] , 
the inclusive probability summed over all such states AD is insensitive to short distances 
and is determined completely by the parameters a, and ry*. The inclusive probabilities 
for producing an atom and a deep dimer from an atom-dimer collision or from the collision 
of three atoms in the state labelled ( J = 0, n) are 

- ■ (50b) 

The factor S21 in Eq- (50a) is the amplitude for an incoming atom-dimer scattering state 
to be transmitted through the scaling region into a hyperradial wave in the scale-invariant 
region. The factor 1/(1 — siie^*^*°~^''*) takes into account an arbitrary number of rcficctions 
of the hyperradial wave from the short-distance region and then from the scaling region. The 
factor 1 — e'^"^* is the probability that a hyperradial wave incident on the short-distance 
region will not be reflected and will therefore ultimately emerge as a scattering state of an 
atom and a deep dimer. The analog of the unitarity condition in Eq. (15) for the case 77* > 
is 

n=3 AD 

If we insert the S-matrix elements in Eqs. (49) and the probability in Eq. (50a), we can use 
the unitarity conditions for the matrix Sij in Eq. (21) to show that Eq. (51) is automatically 
satisfied. 

The hyperangular average of the 3-body recombination rate into the shallow dimer is 
obtained by squaring the S-matrix element in Eq. (49b), summing over the index n labelling 
the 3-atom states, and then multiplying by a kinematic factor. Using the unitarity condition 
in Eq. (51), the expression for the S-matrix element in Eq. (49a), and the expression for the 
probability in Eq. (50a), this can be expressed as 

(l-e-^'?*)|si2(x)p \ ha^ 

(52) 

This reduces to Eq. (22) in the limit r/* 0. The hyperangular average of the inclusive 3- 
atom recombination rate into deep dimers is given by summing the probability in Eq. (50b) 
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over the index n labelling the 3-atom states and then multiplying by a kinematic factor. 
Using the unitarity of the matrix Sij {x) , the rate can be expressed as 



144V37r2(l - e-^^*)(l - |sii(a;)p - 


|si2(a;)P) ha"^ 


x^\l- sii(x)e2'^*o-2'?*|2 


m 



(53) 



This contribution vanishes in the limit 7]^ — > 0. The recombination rates in Eqs. (52) and 
(53) are completely determined by the same universal scaling functions Sii(x), Si2{x), and 
•522(2^) that determine the 3-body recombination rate into shallow dimers in the case 77* = 0. 



B. Analytic results at threshold 



Analytic expressions for the 3-body recombination rates -ft'shaiiow(O) = K^'^\0) and 
-^deep(O) at the threshold E = can be obtained by inserting the analytic results for the 
universal scaling functions su{x), 512(2;), and 522(2^) in Eqs. (34) into Eqs. (52) and (53). 
The 3-body recombination rates at threshold into the shallow dimer and into deep dimers 
are 



7687r2(47r-3v^)(sin2[soln(a/a*o)] +sinh2 77*) x 

-f^ shallow (Uj — —27 : ^? 1 / / ) I04aj 

smh (ttso + V*) + cos^[so ln(a/a*o)J 
. 381.^(1. -3v^),-oth(...,)..nh(2„ ) fto^ 
sinh (ttso + V*) + cos^fso ln(a/a*o)] 

Since sinh^(7rso-l-r;*) > 139 is so large, these expressions can be approximated with errors of 
less than 1% by omitting the cos^ terms in the denominators. To within the same accuracy, 
we can approximate sinh(7rso + ?7*), sinh(7rso), and cosh(7rso) by exponentials to get the 
simpler expressions 

^shaiiow(O) fti eCmaxC"^''* (sin^[soln(a/a*o)] +sinh^77*) ^aV^, (55a) 
i^deep(O) ~ |C„,ax(l-e-^''*)^aVm, (55b) 

where Cmax ~ 67.1 is given in Eq. (32). 



C. Three-body recombination rate 

In this subsection, we use the universal scaling functions calculated in the previous section 
to calculate the 3-body recombination rate into deep dimers and the J — contribution to 

the 3-body recombination rate into the shallow dimer. The effects of deep dimers on the 
contributions to the 3-body recombination rate into the shallow dimer from J > 1 are 
suppressed in the zero-range limit. 

In Eq. (52), the J — contribution K^^\E) to the 3-body recombination rate into the 
shallow dimer is expressed in terms of the universal scaling functions Sii{x), 512(2:), and 
522(2;). In Fig. 6, it is shown as a function of the collision energy E for 6 values of ^*o- 0, 
^TT, ^TT, ^TT, ^TT, and ||7r. The left and right panels are for 77^,0 = 0.1 and 0.2, respectively. 
Comparing with Fig. 3 for 77*0 = 0, we see that the local minimum for 9^q between 
and 0.124 oq becomes less pronounced as 77* increases. When 77* is large enough, the local 
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FIG. 6: (Color online) The J = contribution to the 3-body recombination rate into the shallow 
dimer for cases in which there are one or more deep dimers. The hyperangular average K^^\E) is 
shown as a function of the collision energy E or 6 values of 6**0: ^vr (upper solid line), ^vr and ^vr 
(upper and lower dashed lines), |^7r and ^vr (upper and lower dash-dotted lines), or vr (lower 
solid line). The left and right panels are for 77* = 0.1 and ry* = 0.2, respectively. 



minimum disappears altogether and K^^\E) becomes a monotonically decreasing function 
oiE. 

In Eq. (53), the 3-body recombination rate into deep dimers is expressed in terms of the 
universal scaling functions sii(x) and Si2(a;). In Fig. 7, the maximum and minimum values 
of K^eepiE) with respect to variations of 6**0 are shown for three values of rj^Q-. 0.01, 0.1, and 
0.2. The differences between the maximum and minimum values are only visible at energies 
greater than Ed- The 3-body recombination rate into deep dimers is insensitive to the value 
of 6*^,0, because the dependence on in Eq. (53) is suppressed by a factor of Sii(x), which 
remains small for x < 10. For any value of 77^,, K^eep{E) decreases monotonically with E. If 
there were no relevant length scales at high energy, dimensional analysis would imply that 
Kdeep{E) should be proportional to E^"^ at large E. Our results at the largest values of E 
are compatible with the approach to this simple scaling behavior. 



VII. APPLICATION TO ^^^CS ATOMS 

The Innsbruck group has carried out measurements of the 3-body recombination rate 
for ultracold ^^'^Cs atoms^ in the \f = 3,mj = +3) hyperfine state [16]. By varying the 
magnetic field from to 150 G, they were able to change the scattering length from —2500 ao 
through to +1600 oq. In this range of magnetic field, the \ f = 3,m/ = +3) state is the 
lowest hyperfine state, so 2-body losses are energetically forbidden. Thus, the dominant 



A convenient conversion constant for ^^■^Cs atoms is h^/m = 1.30339 K ag. 
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FIG. 7: (Color online) Three-body recombination rate into deep dimers. The maximum and 
minimum values of Kf^^epiE) with respect to variations of 0*o are shown as functions of the collision 
energy E for three values of 77* : 0.01, 0.1, and 0.2. 



loss mechanism is 3-body recombination. The van der Waals length scale for Cs atoms 
is {mCQ/ti^Y^^ ~ 200ao. The range of scattering lengths studied by the Innsbruck group 
includes a universal region of large negative a and a universal region of large positive a 
separated by a nonuniversal region of small \a\. In the two regions of large scattering length, 
few-body physics should be universal and characterized by 3-body parameters and 77=,, that 
may be different in each universal region. An interesting open question is whether there is 
any relation between the 3-body parameters and for different universal regions. 

In the region of negative a, the Innsbruck group measured the loss rate constant L3 as 
a function of a at three different temperatures: T = 10 nK, 200 nK, and 250 nK. They 
observed a dramatic enhancement of the loss rate for a near —850 oq. At T = 10 nK, the 
loss rate as a function of a can be fitted rather well by the zero-range formula for T = in 
Ref. [15] with parameters a'^ = —850(20) Qq and 77* = 0.06(1), where a'^ denotes the negative 
scattering length for which there is an Efimov resonance at the 3-atom threshold. Thus, the 
large enhancement in the loss rate can be explained by the resonant enhancement from an 
Efimov trimer near the 3-atom threshold. More recently, the Innsbruck group has measured 
the position of the maximum loss rate as a function of the temperature [42] . Its behavior as 
a function of temperature can be explained at least qualitatively by the dependence of the 
binding energy and width of the Efimov resonance on the scattering length [20, 21]. 

In the region of positive a, the Innsbruck group measured the loss rate constant L3 at 
T = 200 nK for values of a ranging from 62.3 ao to 1228 oq. Their results in this region are 
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shown in Fig. 8. The vertical axis is the recombination length ps defined by 



They observed a local minimum in the loss rate for a near 200 a^. This value is near the van 
der Waals length scale {mCQ/h^Y^'^ ~ 200 ao, so range corrections may be large near the 
minimum. The zero-range prediction for the loss rate at T = is a = i^(0)/6, where i^(0) 
is the sum of ii'shaiiow(O) and K^^e^i^) in Eqs. (54). To an approximation of better than 1%, 
they can be approximated by the expressions in Eqs. (55). By fitting the data for a > 500 ao 
to this expression, they obtained a+ = 1060(70) Oq for the scattering length e'^/*^^*°^a*o at 
which the coefficient of achieves its maximum value. Universality would then imply that 
the minimum should be at a*o = 223(15) oq. The fit was insensitive to the value of rj^ and 
yielded only the upper bound r]^ < 0.2. The Innsbruck group also determined the location 
of the minimum directly by measuring the fraction of atoms that were lost after a fixed time. 
The result was amin = 210(10) oq, which is consistent with the value obtained by fitting the 
data for a > 500 Oq. 

We now consider whether our results arc applicable to this system of "^'^'^Cs atoms. 
In the experiments with positive scattering length, the typical peak number density was 
5 X 10^^ cm~^. The corresponding critical temperature for Bose- Einstein condensation is 
Tc 160 nK. Thus, the temperature T — 200 nK was not far above Tc- We ignore this 
complication and calculate thermal averages using the Maxwell-Boltzmann distribution as in 
Eq. (6) instead of the Bose-Einstein distribution. When a = a^nin, the zero-range prediction 
for the binding energy Ejj of the shallow dimer in Eq. (7) gives 3 x 10'^ nK, so T/Ej:, ^ 0.007. 
The binding energy Ejj decreases to about 9 x 10^ nK at a = 1228 ao, so T/ E^ ^ 0.22. Thus 
the temperature T — 200 nK is safely in the region T < 30 Ed in which the thermal average 
in Eq. (6) can be calculated accurately using the scaling functions calculated in Section IV. 

In Fig. 8, we compare the zero-range predictions for the recombination length p3 defined 
in Eq. (56) with the Innsbruck data for T = 200 nK in the a > region. A blow-up of the 
region < a < 400 ao is shown in the left panel of Fig. 9. We assume niost = 3, so that 
L3 = 3a. We set a*o = 210 ao, which is the central vahic of amin measured in the Innsbruck 
experiment. In Fig. 8, we plot pa as a function of a for T = 200 nK and three values of 1]^: 
0, 0.1 and 0.2. The predictions for nonzero 77* are larger than those for r;* = in the region 
a < 500 ao and smaller in the region a > 500 ao- We extend the zero-range predictions 
down to a = 100 ao, which is smaller than the effective range of Cs atoms, even though 
range corrections are expected to be 100% at such small values of a. The predictions for 
rj^ = give a good fit to the data only in the region near 500 ao where the predictions are 
insensitive to The prediction for 77* = 0.2 gives a remarkably good fit to the data for 
a ranging from the highest point at 1228 ao all the way down to 202.7 ao. As a decreases 
below 200 ao, the prediction for r)* — 0.2 continues to decrease monotonically while the data 
seems to have a local maximum near 137.3 ao. The best fit to the data for a in the range 
from 202.7 ao to 1228 ao is 77* = 0.18, which gives a P^r data point of 1.22. 

In the right panel of Fig. 9, we show the zero-range prediction for the lower temperature 
of T = 40 nK. For rj^ — 0, the local minimal of at 210 ao becomes much deeper at 40 nK. 
The predictions for 77^, = 0.1 and 7/^, = 0.2 at 40 nK are essentially indistinguishable from 
the predictions at 200 nK in the region a < 600 ao. 

In Fig. 10, we illustrate the temperature dependence of the recombination length. We 
show p3 for ?7* = 0.2 as a function of a for three temperatures: 0, 40 nK, and 200 nK. 
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FIG. 8: (Color online) The 3-body recombination length p3 for ^^^Cs atoms as a function of a for 
T = 200 nK. The data points are from Ref. [16]. The curves are the zero-range predictions for 
three values of (solid line), 0.1 (dashed line), and 0.2 (dotted line). 




FIG. 9: (Color online) The 3-body recombination length for ^'^^Cs atoms as a function of a 
for T = 200 nK (left panel) and T = 40 nK (right panel). The data points for T = 200 nK are 
from Ref. [16]. The curves are the zero-range predictions for three values of ry*: (solid lines), 0.1 
(dashed lines), and 0.2 (dotted lines). 



Also shown is the Innsbruck data for 200 nK. For a < 500 aO, the difference between the 
predictions at T = and 200 nK is very small. At larger values of a, the prediction is 
significantly lower at 200 nK. The prediction for T = 200 nK fits the data much better 
than the prediction for T = 0. By comparing Figs. 8 and 10, we can see that increasing 
and increasing T both have the effect of decreasing the prediction for ps for a > 500 oq 
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FIG. 10: (Color online) The 3-body recombination length p3 for ^^'^Cs atoms as a function of a for 
?7* = 0.2. The curves are the zero-range predictions for three values of the temperature: T = 
(solid line), 40 nK (dashed line), and 200 nK (dotted line). The data points are for T = 200 nK 
from Ref. [16]. 

without changing the prediction for smaller values of a. Thus accurate measurements of the 
temperature are essential if the data in this region is to be used to determine 7]^,. 

VIII. COMPARISON WITH PREVIOUS WORK 

In this section, we compare our results with previous calculations of the 3-body recom- 
bination rate for identical bosons. The previous results include calculations using specific 
models for the interactions between the atoms as well as calculations based on the univer- 
sality of atoms with large scattering length. 

Esry and D'Incao have calculated the 3-body recombination rate at nonzero temperature 
by solving the 3-body Schrodinger equation for identical bosons interacting through the 
simple potential V{r) = D / cosh^(r/ro) [17, 18]. In Ref. [18], the depth parameter D was 
used to vary the scattering length from +oo to — oo in the region with one 2-body bound 
state and from +oo to — oo in the region with two bound states. This range of D includes 
three universal regions with large scattering length. The range parameter rg determines 
the two 3-body parameters and rj^, for each of those universal regions. To apply this 
simple model to the Innsbruck experiment on ^^^Cs atoms, Esry and DTncao tuned tq to 
fit the position of the recombination resonance at a = —850 Oq. The predicted height of 
this resonance is larger than that measured by the Innsbruck group. Since t]^ is sensitive 
to the height and width of the resonance, tuning tq to get the correct value of gives too 
small a value for r/*. Esry and DTncao also compared the predictions of their model to the 
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Innsbruck data on ^'''^Cs atoms with a large positive scattering length. The model does not 
give a good quantitative fit to the 200 nK data for large positive scattering length. 

Massignan and Stoof have calculated the 3-body recombination rate at nonzero tem- 
perature by solving the Skorniakov-Ter-Martirosian equation in a scattering model for the 
lowest hyperfine spin state of ^^^Cs atoms [19]. Their scattering model is defined by a 
2-body T-matrix element with 4 parameters that models the Feshbach resonance of ^^^Cs 
atoms centered at B = —11 G. The 4 parameters consist of the large background scat- 
tering length Obg = 1800 ao and 3 Feshbach resonance parameters: Bq, A5, and A//. 
Their two-body T-matrix also depends on a background effective range parameter re, which 
provides a "physical" ultraviolet cutoff and guarantees a well-behaved three-body integral 
equation. Massignan and Stoof took to be a function of the Feshbach resonance param- 
eters: Te = — 2?i^/(mabgA/xAi?). The effective range is defined to be twice the coefficient 
of k"^ in the low-momentum expansion of A; cot 5 (/c), where 5{k) is the S-wave phase shift. 
In their model, the effective range is always small and negative, varying from 2re at the 
resonance to —0.3 far from the resonance. An estimate of the true background effec- 
tive range can be obtained by exploiting the fact that the interatomic potential is —Cq/t^ 
at long distances and using Obg and the van der Waals coefficient Cq for Cs atoms as in- 
puts [43, 44]: Tbg ~ +250 oq. Thus the model of Massignan and Stoof may not provide 
an accurate representation of the 2-body problem. For T = 10 nK, the prediction for the 
recombination length p3 at the peak of the resonance in their model is too large by 40%, 
which is about 4 standard deviations. In the nonuniversal region of small scattering length, 
the model predicts a minimum in the 3-body recombination rate near —150 Oq that was not 
observed in the experiment. For positive scattering length, the model does not give a good 
quantitative fit to the 200 nK data. It predicts the local minimum to be at 240 aoi which is a 
little larger than the observed position 210(10) ao. At the minimum, the prediction for pa is 
less than half the measured value, which is too small by more than 10 standard deviations. 
At a = 1228 ao, their prediction for the recombination length p3 is too large by 10%, which 
is about 2 standard deviations. 

Lee, Kohlcr, and Julienne have calculated the 3-body recombination rate at T = using 
a coupled-channel potential model for the relevant hyperfine spin states of ^^"^Cs atoms 
[45]. Their model had 6 fitting parameters: two interchannel coupling parameters that 
were tuned to fit the position and width of the Feshbach resonance and 4 parameters in a 
separable interatomic potential that were tuned to fit the background scattering length, the 
background effective range, and the binding energies of the two shallowest deep dimers [45]. 
They compared their zero-temperature calculations with the measurements of the Innsbruck 
group at nonzero temperature. The agreement between their calculation and the Innsbruck 
data is not very good even at a qualitative level. Their prediction for the magnetic field of 
the recombination resonance in the a < region is 5.1 G compared with the measured value 
of 7.4 G [42]. Away from the resonance, their predictions for (ps)^ are sometimes lower than 
the 10 nK data by about 2 orders of magnitude. In the a > region, their predictions for 
{pzY were below the 200 nK data by at least one standard deviation at all points and by 
many standard deviations at some of the points. Their prediction for the position of the 
recombination minimum is about 130 ao, which is much smaller than the observed value 
210(10) ao. 

Lcc, Kohler, and Julienne have criticized universal approaches to the 3-body recombina- 
tion as being "incomplete" and they claimed that these approaches "could therefore lead to 
unreliable conclusions" [45]. These criticisms were based on misconceptions. In the case of 
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^^^Cs atoms, the two deep dimers that are closest to the threshold have an avoided crossing 
near the Feshbach resonance at —11 G. Lee, Kohler, and Julienne pointed out that as a 
consequence the probabilities Z_i{B) and Z_2{B) for these deep dimers to consist of atoms 
in the lowest hyperfine spin state change dramatically with the magnetic field. They con- 
cluded incorrectly that ii'deep(O) must depend on B not only through the scattering length 
a{B) but also through the strong 5-dependence of the properties of the deep dimers. It 
should be noted that no quantitative evidence for this conclusion was given in Ref. [45]. 
Their conclusion would certainly be true of the exclusive recombination rates into each of 
the deep dimers. What the authors of Ref. [45] did not appreciate is that the inclusive 
recombination rate summed over all deep dimers is much less sensitive to details involving 
short distances, such as level crossings between the deep dimers. The inclusive effects of 
the deep dimers can be taken into account through the 3-body parameter r/*. Universal 
approaches are certainly "incomplete" in the sense that, although they can describe the 
inclusive recombination rate into deep dimers and the exclusive recombination rate into the 
shallow dimer in the case a > 0, they cannot describe the exclusive recombination rates into 
individual deep dimers. However, contrary to the suggestion of Ref. [45], there is no reason 
to expect universal approaches to lead to "unreliable conclusions" for those recombination 
rates that they can describe. 

In Ref. [45], Lee, Kohler, and Julienne referred to universal approaches to the 3-body 
problem as "universal fitting procedures". Assuming that the scattering length is known, 
the zero-range results require fitting the 3-body parameters and r/*. This approximation is 
very economical in the number of fitting parameters. One price that is paid for this economy 
is that the 3-body parameters depend sensitively on the details of the 2-body potential. They 
may be so sensitive that the only practical way to determine them is through experimental 
measurements of 3-body observables. The model used in Ref. [45] has 6 parameters that were 
used to fit results from a more accurate representation of the 2-body problem [46]. Three 
of the parameters are sufficient to reproduce the scattering length a{B) as a function of the 
magnetic field. The other 3 were used to fit the background effective range and the binding 
energies of the two shallowest deep dimers. The authors expected this to be sufficient to 
calculate the 3-body recombination rate accurately. The large errors in their predictions for 
the positions of the recombination resonance for a < and the recombination minimum for 
a > shows that this is not the case. It implies that their model does not correctly predict 
the physical values of /t* for these two universal regions. One way to proceed would be to 
add more and more parameters to the interatomic potential and use them to fit additional 
2-body observables. However and 77* might be too sensitive to the details of interactions 
at short distances to make this approach practical. An alternative that is motivated by the 
universal approach is to add 2 parameters to the interatomic potential and use them to tune 
appropriate 3-body observables in the universal region to their measured values, thus tuning 

and r)^ close to their physical values. 

In Refs. [22-24], the 3-body recombination rate for the case a > has been calcu- 
lated using simplifying assumptions to neglect some of the universal scaling functions. The 
starting point in Ref. [22] was an approximation in which only the first term in the expres- 
sion for K^'^^E) in Eq. (14), which corresponds to the 3-body channel labelled by n = 3, 
was retained. Inserting the zero-range expression for the S-matrix element S^j^^J^aI-^) 
Eq. (19b), this approximation becomes 
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This approximation could have been avoided by using Eq. (22) for K^^\E) as a starting 
point instead. The first simphfying assumption of Ref. [22] was to neglect Sii(x) in the 
denominator. According to the analytic result in Eq. (34a), this function is very small at 
the threshold: |sii(0)| a; 0.002. The assumption of Ref. [22] was that |sii(a;)| remains small 
for larger values of x. This assumption has now been verified by our explicit calculation of 
Sii(a;), which shows that |sii(a;)| < 0.035 for x < 10. With this simplifying assumption, the 
expression for the recombination rate in Eq. (57) can be reduced to 

i^W(E) ^ 1^1^ \s,,{x)e-^'*^ + s,,{x)s,,{x)e^''^\' — . (58) 

The amplitude inside the absolute value sign can be expressed as a linear combination of 
cos ^-nd sin d^o with coefficients that are complex functions of x. The coefficient of sin 
can be made real-valued by multiplying the amplitude by an x-dependent phase that does 
not affect the rate. Thus the approximation in Eq. (58) depends on 3 independent real- valued 
scaling functions. The second simplifying assumption of Ref. [22] was that the coefficients 
of cos^^*o ^-nd sin^*o were relatively real. This assumption was motivated by the existence 
of a deep minimum in K^^\E) at = 2.1 mK in the 3-body recombination rate for "^He 
atoms calculated in Ref. [38]. This behavior can be most easily understood if the square of 
the complex amplitude in Eq. (58) can be approximated by the square of a real amplitude 
that vanishes at = 2.1 mK. Given the second simplifying assumption, the approximation 
in Eq. (58) depends on only two independent scaling functions. In Ref. [22], these two 
scaling functions were constrained by the requirement that they reproduce the result of 
Ref. [38] for K^^\E) for "^He atoms interacting through the HFD-B3-FCI1 potential, which 
corresponds to a — 1.15 a*o- By assuming the absence of a near cancellation between the 
two scaling functions for a — 1.15 a*o, they obtained predictions for the recombination rate 
with surprisingly narrow error bands in the region of a near a*o- In Ref. [23], Shepard 
reduced the error bands to a line by using K^^^E) for ^He atoms interacting through a 
second potential to determine the two scaling functions separately. In Ref. [24] , Platter and 
Shepard removed the second simphfying assumption of Ref. [22] by using K^^\E) for four 
different ^He potentials to determine the three independent real-valued scaling functions in 
Eq. (58). Their approximation to the zero-range result for K^^\E) is accurate to better 
than 5%. 

In Refs. [22] and [24] , the universal scahng functions that were extracted from calculations 
of 3-body recombination for ^Hc atoms were applied to ^^^Cs atoms with a > 0. The effects 
of deep dimers on the 3-body recombination rate into the shallow dimer were obtained by 
making the substitution in Eq. (48) for loga*o in the expression for K^^\E) in Eq. (58). 
Since their approach was unable to determine the function Si2{x) that appears in Eq. (53), 
they approximated the 3-body recombination rate ii'deep(-E') into the deep dimer by its 
value at E = 0, which is given in Eq. (55b). In Ref. [22], the error bands associated with 
undetermined universal scaling functions were small for a < 300 Oq but they widened rapidly 
as a increased, so they were unable to obtain a quantitative fit to the Innsbruck data. In 
Ref. [24], the authors obtained a reasonably good fit to the Innsbruck data for a > 200 oq 
with values of ry* in the range from to 0.01, although their results tended to lie above the 
highest data point at 1228 ao- Onr definitive zero-range results give a significantly better 
fit. There is a qualitative difference in the dependence on 77* between our definitive zero- 
range results and the approximation in Ref. [24]. As rj^ increases, the definitive zero-range 
prediction in the region a > 600 Oq decreases as shown in Fig. 8, while the approximation 
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in Ref. [24] increases. The difference comes primarily from the approximation of Kdccp{E) 
by its threshold value in Ref. [24]. As illustrated in Fig. 7, this is a good approximation for 
£^ < 0.1 Ed, but Kdeep(-E) decreases rapidly with E for energies above Ed. The calculation 
of the energy dependence of K^eepiE) was essential for obtaining the excellent fit to the 
recombination length of ^^^Cs atoms with 77* = 0.2 that is illustrated in Fig. 8 and the left 
panel of Fig. 9. 

IX. SUMMARY AND OUTLOOK 

In this work, we have calculated the universal scaling functions that determine the 3-body 
recombination rate for identical bosons as a function of the collision energy in the zero- 
range limit^. We used the STM equation to calculate the phase shifts for elastic atom-dimer 
scattering above the breakup threshold for angular momenta up to J = 6. For J > 1, the 
contribution K^"^^ (E) to the 3-body recombination rate into the shallow dimer is expressed 
in terms of a real- valued universal scaling function fj{x) in Eq. (17). We determined fj{x) 
for each J > 1 by fitting the imaginary part of the atom-dimer phase shift. For J — 0, the 
contribution K^'^\E) to the 3-body recombination rate into the shallow dimer is expressed 
in terms of three complex- valued universal scaling functions Sii{x), Si2{x), and S22{x) in 
Eq. (22). We determined these scaling functions by fitting the S-wave atom-dimer phase 
shift to Efimov's radial law. The parameters in the zero-range predictions for the 3-body 
recombination are the scattering length and the 3-body parameter a^o- We compared our 
results to previously published calculations of the 3-body recombination rate of "^Hc atoms 
for angular momenta up to J = 3 [38]. We used the 3-body recombination rate at threshold 
calculated in Ref. [38] to determine a*o for ^He atoms. We found good quantitative agreement 
for the 3-body recombination rate K'^'^^E) as a function of energy in the J — 0, 1, and 2 
channels. 

Our results can also be applied to atoms with large scattering length that have deep 
dimers. The deep dimers provide additional 3-body recombination channels and they also 
have an effect on the 3-body recombination into the shallow dimer. An additional 3-body 
parameter 77* is necessary to take into account the inclusive effects of the deep dimers. Re- 
markably, the 3-body recombination rate is completely determined by the same universal 
scaling functions that apply at ?]* = 0. The J = contribution K^^\E) to the 3-body 
recombination rate into the shallow dimer is expressed in terms of the three complex- valued 
universal scaling functions Sii(a;), 512(2;), and S22{x) in Eq. (52). The J > 1 contributions 
K^^\E) are unaffected by the deep dimers and are still given by Eq. (17). The inclusive 
rate K^^epiE) for 3-body recombination rate into deep dimers is expressed in terms of the 
universal scaling functions Sii{x) and 512(2;) in Eq. (53). We compared our results to mea- 
surements by the Innsbruck group of the 3-body recombination rate for ^^^Cs atoms with a 
large positive scattering length at a temperature of 200 nK. We set the 3-body parameter 
a*o equal to the position of the recombination minimum measured by the Innsbruck group 
[42]: a*o = 210 gq. We then used the 3-body parameter 77* to fit the dependence of the 
recombination length on a. We find that 77* — 0.2 gives excellent agreement with the exper- 
imental results for T = 200 nK in the region a > 200 oq. We also found that lowering the 
temperature would have an observable effect only for a > 500 Qq. 



^ Numerical results for the universal scaling functions can be obtained by contacting any of the authors. 
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Our excellent fit to the T = 200 nK data is particular remarkable considering that it was 
obtained by using only the scattering length and two adjustable 3-body parameters. The 
accuracy of this fit presents a challenge to approaches that use more fitting parameters but 
tune them to fit 2-body observables in addition to the scattering length, as in Ref. [45]. 
The 3-body parameters and r]^ associated with a given universal region arc complicated 
functions of the 2-body parameters. The 3-body parameters and rj^ may be too sensitive 
to details at short distances to be determined accurately by fitting a small number of 2-body 
observables. It is possible that the only practical way to determine them accurately is by 
tuning them to fit measurements of 3-body observables. 

Our zero-range results for the 3-body recombination rate have a fractional theoretical 
uncertainty rs/\a\, where Vg is the effective range. Higher accuracy could be obtained by 
calculating range corrections. It is known how to account systematically for these corrections 
using an effective field theory for short-range interactions [11, 13]. A first calculation of 
the range corrections to the recombination rate into the shallow dimer at threshold and 
an application to ^He atoms has already been carried out [41]. In the future, the range 
corrections should be calculated for the recombination rate as a function of the energy. The 
range corrections to the recombination rate into deep dimers should also be calculated. The 
range corrections are particularly important near a recombination minimum. In particular, 
they can shift the position of the minimum [41]. In the case of Cs atoms, the importance of 
range corrections near the recombination minimum is further increased by the fact that the 
minimum is at a scattering length 210 Oq that is not large compared to the effective range. 

Our calculations and those in Refs. [18, 19] have been carried out in the Boltzmann 
region of temperature and density. It would be useful to have zero-range predictions for 
lower temperatures where the Bose-Einstcin distribution has to be used. In the Boltzmann 
region, we only need the hyperangular average K[E) of the 3-body recombination rate, 
which can be determined from the atom-dimer phase shifts using unitarity. In the Bose- 
Einstein region, it is necessary to calculate the 3-body recombination rate directly. This 
requires the calculation of scattering rates with 3 atoms in the initial state. The relevant 
S-matrix elements are expressed in terms of universal scaling functions ,sii(a;), Si2(a;), Si„(a;), 
and S2n{x) in Eq. (19b). The scaling functions Sin{x) and S2n{x) would have to be calculated 
for those 3-atom states labelled by n = 3, 4, . . . that are most important near the 3-atom 
threshold. 

It is also important to calculate the definitive zero-range predictions for the 3-body re- 
combination rate for negative scattering length. In this case, recombination can occur only 
into deep dimers. Previous calculations based on the universal approach have used the adi- 
abatic hyperspherical approximation [20] or a resonance approximation [21]. The relevant 
universal scaling functions can be determined by calculating 3-atom elastic scattering rates 
near the 3-atom threshold. 
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